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ABSTRACT 

The role of matter enhanced resonant neutrino oscillations in reviving a stalled 
shock in a type-II supernova through delayed neutrino heating is investigated. The 
extent of neutrino heating is estimated for the allowed possibility of complete flavour 
conversion self-consistently with the changes in nuclear equilibrium. The average in¬ 
ternal energy per nucleon is substantially increased indicating the possiblty of a robust 
explosion. 
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The theory of supernova explosions has seen a lot of progress in the last hfteen 
years but it is still not able to deliver an explosion with the right energy. Though most 
of the ingredients developed for the scenario of gravitational collapse and shock propa¬ 
gation after core rebounding are believed to be correct and some have got conhrmation 
through the detection of neutrinos of SN1987A, the strong dissipation of shock energy 
through nuclear dissociation and neutrino losses make the shock either fail to reach 
the edge of the core or emerge, after neutrino heating, with energies less than 10®^ 
ergs, a factor of 3 to 4 less than observed values. Different mechanisms like improved 
pre-supernova conditions, soft equation of state, general relativity at high densities, 
convection and improved neutrino physics have all been invoked to solve this problem. 
In this letter we describe a novel scenario by coupling neutrino flavour conversion in 
matter (the MSW mechanism) with the delayed neutrino heating that can deliver ex¬ 
tra energy to matter behind the shock to generate a more robust explosion. Fuller 
et. al. m noted this possibility of conversion of mu/tau neutrinos to electron type 
between the neutrino-sphere and the stalled shock for some ranges of neutrino mixing 
parameters and estimated the extra heating and the explosion energy. We incorporate 
this heating to a semianalytic evolutionary calculation of the thermodynamic condi¬ 
tions behind the shock front of a 25Mq star and demonstrate explicitly the change 
in nuclear statistical equilibrium (NSE). We also calculate the gain in internal energy 
achieved by the increased heating induced by neutrino flavour oscillations. The results 
are encouraging. 

At present there are two possible scenarios for supernova explosion. For stars in 
the mass range 8Mq < M < 15Mq, under some very special conditions on the size 
and structure of the core and the equation of state, the shock continues its outward 
propagation and the star explodes within some tens of milliseconds after the beginning 
of collapse. This is the prompt explosion scenario 0. For more masssive stars the 
energy of the shock gets dissipated in dissociating nuclei and producing uu pairs and 
the shock stalls at a radius of a few hundred kms and becomes an accretion shock. It 
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is subsequently revived by the heating caused by neutrinos from the neutrino-sphere. 
This is the delayed explosion mechanism P, This late time neutrino heating behind 
the shock is caused by their absorption reactions on the nuclei as well as free nucleons 
and by charged and neutral current scattering reactions. Ray and Kar [^] considered 
a schematic model for delayed neutrino heating with matter in the form of a typical 
heavy nucleus ^^Fe, a-particles, neutrons and protons. They calculated the amount of 
neutrino heating treating changes in the NSE behind the shock front in a self-consistent 
fashion. Subsequently, Haxton [Q pointed out the necessity of including the neutral 
current inelastic scattering of nuclei to giant resonance states in calculating the neutrino 
heating rates. Although all three flavours participate in this, the and Uj. neutrinos 
contribute more to the heating by the above process because of their higher equilibrium 
temperature. In |P, the charged current cross-sections of z/gS and z^s on nuclei for some 
typical neutrino temperatures are also tabulated. Ray and Kar [0] incorporated these 
cross-sections in their semianalytic model and computed the thermodynamic state 
of of matter behind the shock as a function of time and estimated the heating rates. 
Fuller et. al. pointed out that matter enhanced oscillations of neutrinos in the 
region between the neutrino-sphere and the stalled shock can increase the heating rate 
appreciably and may result in a delayed explosion with an energy > 10^^ ergs, for a 
cosmologically signihcant or Ur mass of 10 -100 eV and small vacuum mixing angles. 
Due to neutrino flavour mixing z/^s or (which carry a higher average energy) get 
converted to z/gS, the upshot being the production of higher energy z/gS which can heat 
the shock more effectively. 

We re-examine the situation modelled in |0 with the addition of neutrino heating 
in the presence of matter induced neutrino oscillations. We consider a thin layer of 
matter behind the shock, at a radius Rm, which gets heated by the neutrinos coming 
from the neutrino-sphere situated at a distance Ry from the center. Matter enhanced 
resonant neutrino oscillations have been considered extensively in the context of the 
solar neutrino problem IQ. The intrinsic neutrino properties needed for the purpose are 
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neutrino mass squared difference, Am^, and the mixing angle in vacuum, 6v- Oscilla¬ 
tion of neutrinos in matter is different from that in vacuum because interactions modify 
their dispersion relation and consequently they develop an effective mass dependent 
on the matter density. The mixing angle in matter can be expressed as 


tan29M = taii29v/{1 — E^/Ea) 


( 1 ) 


El, is the neutrino energy; Ea is given by Ea = cos29v/2\/2GFnei where Ue 

denotes the electron density. As is seen from (2) tan26'M goes through a resonance 
when Ey = Ea- From this one gets an expression for the resonance density as |I|, 


pres = (2.108 X 10^°gm/cc) ^ 


0.5 \ f Am'^ cos29v\ /lMeV\ 


V 1600eV2 J \ El, 


■) 


( 2 ) 


At resonance, 9m is vr/d and one gets maximal flavour mixing. The effectiveness of 
neutrino oscillations in increasing the heating rate depends crucially on whether a 
resonance is encountered between the neutrino-sphere and the shock front or not. If 
a neutrino of flavour / encounters a resonance between the neutrino-sphere and the 
shock then its probabiity to remain a Uf at the position of the shock is given as. 


P{uf Uf) = 0.5 -F (0.5 - Pj) cos 29r^ cos 29r^ (3) 


where 9r^ {(^Rm) is the neutrino mixing angle in matter at the position of the neutrino- 
sphere (shock). Pj is the non-adiabatic transition probability between the two neutrino 
states and can be expressed in the Landau-Zener approximation as, P j = cxp^—Eria/E,,) 
where. 


, , sin^26'y (Am?\ f rres\ 

£„. = (0,026MeV)^ ^ (ifS) 


(4) 


Tres is the resouauce radius. Eria can be obtained using (||) and a density prohle 0 
p=(10^^gm/cc)(r/lcm)“^. The above equations are for oscillation between two neu¬ 
trino flavours. In what follows we take these two flavours as the electron and the tan 
neutrino. A neutrino passing through the resonance density can undergo complete 
flavour transformation for suitable values of the parameters. A total conversion of the 
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which carry a higher energy, to Ve between the neutrino-sphere and the shock will 
generate maximum heating. Ideally, one should couple the neutrino flavour transforma¬ 
tions with the hydrodynamics of the shock and calculate the probabilities of conversion 
between the neutrino-sphere and the shock at each instant and And the altered heat¬ 
ing rates and the consequent changes in NSE self-consistently. In this letter we have 
estimated the heating rates for full flavour conversion between the neutrino-sphere and 
the shock. 

We consider matter consisting of nucleons as well as nuclei and incorporate charge 
current cross-sections on nuclei from [Q, where these are listed upto a temperature of 
6 MeV. If flavour conversion takes place, the z/g temperature can be higher. A simple- 
minded extrapolation yields a much higher cross-section and consequently a higher 
heating. We take a conservative standpoint and use the cross-section given for 6 MeV as 
the value at higher energies. We also incorporate the neutral current inelastic scattering 
on nuclei. As pointed out in [§] and later shown explicitly in [^, such reactions can 
increase the shock heating rates substantially. However, these processes being flavour 
independent, the heating due to them is unaffected by neutrino oscillations. The energy 
absorbed by matter behind the shock front/gm/sec is given by, E = Ei — E 2 + Egcatt, 
El is the heating rate due to charged current absorption of z/, z/ on free nucleons and 
nuclei. E 2 is the rate of energy loss due to neutrino radiation at matter temperature 
Tm and Escatt is due to u-e scattering and is T^/Ty times the difference of the first two 
terms @] plus the contribution due to inelastic neutral current scattering on nuclei 0 . 
With complete flavour transformation these terms are different from what they were in 
the absence of oscillation ii , Ty^ in these terms being replaced by Ty^. The resulting 
increase in E is quite significant. The energy absorbed first goes into dissociating the 
iron nuclei into a-particles via ^^Ee ^ 13a -|- 4rz,. The a-particles then break up into 
protons and neutrons through a ^ 2p + 2n. The Q-values per nucleon for these are 
2.2 MeV and 7.075 MeV respectively. The exact nuclear composition at each instant 
of time is governed by the Saha equations 0. The change in entropy AS due to the 
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absorption of energy and changes in nuclear composition in time At is given, in the 
limit of zero mass accretion, by 

TmAS = EAt - 2.2AXFe - 7.075AX„ (5) 

In this Xi denotes the relative fraction of species i, AXi is the change in time At. For 
non-zero mass accretion there will be additional terms in (|^) due to the kinetic energy 
of the accreting matter plus the terms due to compositional changes of the matter as it 
moves across the shock front The corresponding change in temperature AT is obtained 
as in 1^. The internal energy of the matter per nucleon tint has contributions coming 
from baryons, electrons and photons. The success or failure of this phase of heating 
depends on the comparative values of eint and = GMc/Rc, where Rc denotes the 
bifurcation radius and Me is the mass included within R^. 

For illustrative purposes numerical calculations are performed for a 25Mq star p] 
with an initial iron core mass of 1.37M0. The electron neutrino luminosities for 
which calculations are performed are 4 x 10^^ and 2 x 10^^ ergs/sec with the corre¬ 
sponding temperatures of neutrino-spheres as 4.45 and 5.3 MeV. The luminosities of 
neutrinos and antineutrinos of all flavours are assumed to be the same. The z/g-sphere 
radius, Ry, is 30 km. The temperatures of the fi and r neutrino-spheres are taken as 
10 MeV. The shock stalls initially at a radius of 480 kms. It then recedes upto 189 
kms and moves forward again |0]. We consider non-zero mass accretion with the rate 
0.1 Mq/sbc. 

Demanding that a resonance is achieved between the neutrino-sphere and the shock 
so that Ry < Tres < Rm, an idea regarding the bounds on Am^ can be obtained under 
the assumption of small mixing angles. For complete flavour conversion one has to 
ensure that not only is a resonance attained, but P{i^f —>■ Vf) as given by is zero. 
(99 - 100)% flavour conversion of all neutrinos in the energy range 1-100 MeV requires 
1.13 X lO^eV^ < Am^ < 1.62 x 10%V^, 1.83 x 10“^ < sin%0y/cos20y < 7.85 x lO'T 
This is the most stringent bound on the parameter values. However, heating is most 
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effective for neutrino energies 10-35 MeV. For such neutrinos the allowed area in the 
parameter space is 2.95 x lO^eV'^ < Am? < lO^eld^, 1.45 x 10“® < sin^26'y/cos 26*y < 
3.39 X 10“^. This is consistent with the cosmological limit of < 100 eV 


The nonadiabatic MSW solution to the solar neutrino problem, which requires Am? ~ 
6 X 10~^eV^ and sin^ 29v ~ 7 x 10“^ is due to Ue — oscillations in this framework. 

Figure la shows the time evolution of the fraction of a-particles, neutrons and 
protons with and without oscillation for = 4 X 1052 ergs/sec. The starting point of 
our calculations t=0.417 sec is the time at which the shock starts retreating, where t=0 
corresponds to core bounce. By this time iron is completely dissociated for the 25Mq 
star considered. The energy delivered to the matter goes to dissociate the a-particles 
present. As is seen from hgure la, in the presence of neutrino flavour conversion, 
dissociation of a’s takes place at a faster rate because of the enhanced heating. This 
further raises the heating rate due to an increased free nucleon fraction. Figure lb 
shows the total internal energy per nucleon as a function of time for two representative 
luminosities = 4 x 10®^ ergs/sec, 2 x 10^^ ergs/sec. For a 25Mq star, for the two 
chosen luminosities, even without oscillation eint is greater than the critical energy Ccr- 
For Ly^ = 4 X 10®^ ergs/sec, neutrino oscillation increases it by 61.33% at the end of 
this heating phase. In case of Ly^ = 2x 1052 ergs/sec, was just above the critical 
value but complete flavour conversion raises it by 56.37%. 

r 2 [ indicate that this phase of shock revival is followed by 


The results of or 

the creation of a quasi-vacuum, dominated by radiation, which ends with an eventual 
expulsion of matter. Our initial results indicate that resonant neutrino oscillations may 
indeed be the important ingredient, hitherto missing, in the delayed neutrino heating 
scenario to give a stronger shock in conformity with observation. One needs to extend 
these calculations to later times as well as incorporate the physics of neutrino oscil¬ 
lations in hydrodynamic codes. We are also studying the effect of these oscillations 


on a 18Mq star |]I^ with neutrino luminosities as given in [|I^] for its importance in 
connection to SN1987A. Earlier calculations for a I 8 M 0 star without neutrino oscilla- 
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tions failed to strengthen the shock adequately [0 . Our preliminary calculations show 
that with matter induced neutrino oscillation the internal energy becomes greater than 
the critical energy in this case as well. In conclusion, we stress that resonant neutrino 
oscillations open up exciting possibilities for type-II supernove through the delayed 
neutrino heating mechanism and merit detailed examination. 
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FIGURE CAPTION 


Fig.l. Time evolution of (a). Xp, X„, X^ for a 25Mq star with = 30 km. The solid 
(dashed) curves are with (without) neutrino oscillations, (b). The internal energy per 
nucleon e^ni of a mass shell behind the shock front. Curves I (II) are for luminosities 
= 2 X 10^^ {Li, = 4 X 10^^) ergs/sec. Also shown is the critical energy ecr- 
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